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ABSTRACT 

Exact solutions of conformal turbulence restricted on a upper half plane are 
obtained. We show that the inertial range of homogeneous and isotropic turbulence 
with constant enstrophy flux develops in a distant region from the boundary. Thus 
in the presence of an anisotropic boundary, these exact solutions of turbulence 
generalize Kolmogorov's solution consistently and differ from the Polyakov's bulk 
case which requires a fine tunning of coefficients. The simplest solution in our case 
is given by the minimal model of p = 2, g = 33 and moreover we find a fixed point 
of solutions when p, q become large. 
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In the theory of turbulence, Kolmogorov's theory of the inertial range provides 
a simple but powerful tool in obtaining the scaling laws of locally isotropic tur- 
bulence!^' Nevertheless, the microscopic, field theoretic foundation of the scaling 
law so far has been lacking. Recently, Polyakov has proposed a completely new 
approach to this problem based on the assumption of conformal invariance of ve- 
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locity correlators in the case of two dimensional fully developed turbulence. ' He 
took the non-unitary minimal conformal field theory as a field theoretic framework 
for turbulence and obtained exact solutions of turbulence with specific scaling be- 
haviors, different from the Kolmogorov's one, which seem to fit better with exper- 
imental results!*' However, due to the nature of non-unitarity, one point functions 
of operators in the theory do not necessarily vanish unlike the unitary case. The 
arbitrariness of one point functions in the Polyakov's approach makes the theory 
largely undetermined and remains as an open problem. 

In this letter, we show that one point functions are determined by boundary 
conditions, as suggested by Polyakov, in particular when they are confined on a up- 
per half plane. This allows us to give a field theoretic account for the Kolmogorov's 

[51 

hypothesis of inertial range. We show that the inertial range of homogeneous and 
isotropic turbulence with constant enstrophy flux develops in a distant region from 
the boundary. Moreover we argue that exact solutions obtained by Polyakov and 
others'" for the bulk case require an unnatural flne tunning of coefficients. With- 
out such a fine tunning, we obtained series of solutions of conformal turbulence 
which are different from the bulk case. The simplest solution in our case is given 
by the minimal model of p = 2, g = 33. More interestingly, we find a fixed point of 
conformal solutions when p, q become very large whose energy spectrum exponent 
is -67/15 ~ -4.4667. 

The basic assumption of Polyakov's conformal turbulence is to identify the 
stream function ip, arising from the velocity vector Vi of incompressible fluid such 
that Vi — eijdjip, with a certain primary operator of a non-unitary minimal model 
which, in the inviscid and static case, satisfies the Navier-Stokes equation as well 
as the constant enstrophy condition. The use of two dimensional conformal field 
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thcory(2-d CFT) allows us to compute exactly the velocity correlation functions, 
thereby the kinetic energy of fluid motion. In particular, the velocity two point 
function can be obtained directly from the stream two point function by differen- 
tiation, for example, {vx{zi)vx{z2)) = dy^dy2{ip{zi)i>{z2)); z ^ x + iy. 

In principle, the two point function of stream function can be obtained exactly. 
However, in practice, for the brevity of our discussion and also in order to avoid 

the difficulty in solving for exact solutions, we consider only the short distance 
behavior of the velocity two point function which can be obtained easily from the 
operator product expansion (OPE) of two V^'s; 

(^(.i)^(.2)) = E ^fe , _,%-2A,^ + descendents, (1) 

where C^^ are structure constants and (pi are primary operators of conformal 
dimension A^.. The one point functions except for that of the identity operator 
vanish identically when the vacuum is invariant under SL{2,C). However, in 
the presence of large scale structures of scale L, so that the vacuum is no longer 
invariant under SL{2, C), the one point function behaves as 

{(f>i) ^ CiL-^^^i . (2) 

For the unitary case, A^. > and approaches zero when L becomes large. 
However this is no longer true for the non-unitary model where A^. could be 
negative. Moreover non zero values of one point functions modify significantly 
other correlation functions so that they become essential ingredients in defining a 
non-unitary field theory. 

In order to understand the role of the one point function and its physical im- 
plications more clearly, we first note that they can be determined exactly by a 
specific boundary condition.'"' Consider, for example, the two dimensional turbu- 
lence restricted on a upper half plane; Im 2; = y > 0. A typical boundary condition 



3 



is such that the normal component of the velocity vector vanishes at the boundary. 
In order to maintain the spirit of 2-d CFT approach to turbulence on a upper half 
plane, we assume that the conformal symmetry of correlation functions persists in 
the presence of boundary, but now in a reduced form which preserves the bound- 
ary condition. This is tantamount to reducing SL{2, C) to SL{2, R) in the case of 
the global conformal transformation so that only real translations and dilatations 
are involved. The domain under consideration, the upper half plane, is invariant 
under global real conformal transformations of SL{2,R)/Z2 with an added point 
at infinity. This boundary condition also requires correlation functions to be in- 
variant under SL{2,R). In the case of two-point function, one could obtain this 
by forming an SL{2, R) invariant u 

^ ^ {zi - zi){z2 - Z2) 

(Zl - Z2){Z2 - ZlY 

where two points zi and Z2 are given in the upper half plane and zi and Z2 are 
their images in the lower half plane. In general, it is possible to show that for the 
degenerate conformal field theories, the n-point function {(j){z\,'zi) ■ ■ ■ 4>{zn,Zn))b 
can be obtained systematically from the bulk 2n-point function, where (• • de- 
notes correlation functions in the upper half plane!^°' This is so because the n-point 
function {(l){zi,zi) ■ ■ ■ (f){zn, Zn))b satisfies the same differential equation as does the 
bulk 2n point function consisting of charges in the upper half plane as well as their 
images in the lower half plane. In particular, one point functions on a upper half 
plane can be obtained exactly from the known bulk two point functions such that 

mz))b = d^,y-^^^^; y = lmz, (4) 

where rf^^ are arbitrary constants. These constants d^^, as well as the negative 
energy states (A^^ < 0) in non-unitary minimal models, seem to be unphysical at 
first sight. However, they acquire physical meaning when the non-unitary model 
is applied for an effective description of an open physical system. In our case, 
d(j). obtain physical meaning through the velocity correlators and become physical 
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parameters of turbulence. From Eq.(4), for example, we have an average velocity 
profile; 

{vy)b = -dMh = 0; = dy{'iP)i, = -2d^^^y-''^^-\ (5) 

where is a parameter which controls the magnitude of average horizontal veloc- 
ity. Another example where o?^. manifest physical meaning is the kinetic energy of 
turbulent fluid. It is commonly expressed in terms of the energy density at point 
(x, y) given by 

\{vl{x,y))b = j dkE^^ y){k). (6) 

In the momentum space, the energy spectrum E^^^y^^{k) is given in terms of velocity 
correlators: 

E{x,y){k) = ^ y d'x'e'^-^----'-\vfi{x', y')vp{x, y)),. (7) 
For a large k, E(^^ y-^{k) can be obtained directly from Eqs.(4) & (1) such that 

Eix,y){k) - Yl^%k^^''-^^'"^\<l>i)b + descendents 

' . (8) 
= J^'^ti'^^^^^'^^'^'^^^'t'^y'^^^^ + descendents. 

i 

Thus we have obtained an explicit expression of the energy spectrum as a function 
of the distance from the boundary. Unlike the bulk case, the spectrum is made of 
a sum of different powers controlled by A^^ and (i^^. The appearance of too many 
parameters, in general, is not desirable in the statistical description of turbulence. 
However, if we move far away from the boundary, i.e. when y becomes very large, 
only the first few terms in the series of Eq.(8) dominate and the theory is controlled 
by only few parameters. As we will discuss below, this is precisely in agreement 
with the traditional Kolmogorov's idea of the inertial range. 
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The main idea of Kolmogorov is first to assume that with sufficiently large 
Reynolds number, turbulence becomes locally isotropic - i.e. small scale velocity 
field fiuctuations are statistically homogeneous, isotropic, and stationary - regard- 
less of the form of a finite space-time region bounding turbulence. Then Kol- 
mogorov suggested that for locally isotropic turbulence, the energy spectrum in 
the intermediate scale (inertial range) depends only on the constant flux of energy 

[13] 

(or enstrophy for the two dimensional turbulence ). In our case, we note that 
for sufficiently large y, conformal turbulence becomes locally isotropic. This is so 
because even in the presence of the anisotropic half plane boundary, the small scale 
fluctuations of velocity fleld become isotropic when y becomes large. In order to 
see this explicitly, one could compare for example a small scale velocity fiuctuations 
{vx{x,y)vx{x,y + e))i) with its 90° rotation {vy{x,y)vy{x + e,y))ii;e <C 1. In the 
leading order, they both agree to 



2CJ,(2A^- 

{vj;{x, y)vx{x, y + e))^ ^ {vy{x, y)vy{x + e, y))b ~ oA^.4A^-2A^+2 " ' (9) 

where (j) is the minimal dimension operator of conformal dimension coming 
from the operator product On the other hand, the constant enstrophy 

condition which follows from the r-independence of {w{r)w{0)) , can be met by a 
single condition when y is large 

+ - = -3. (10) 

where x is the minimal dimenson operator coming from the OPE of [0] * In case 
— 0, the constant enstrophy condition changes into 

A,^ + A^ = -3. (11) 

Therefore, we contend that the inertial range forms in a distant region from the 



* The minimal dimensional operator x from the OPE of [V'][</>] is not necessarily the same 
as the one coming from [-0] [■</>] [f/^'J , which were falsely identified in Ref.[3]. In many cases, 
they are equal but as in the Table 1, M^^y^ M^^n are exceptions. In the case of Ms^ya, ij} = 
^(2,26), 4> = *(i,i5), X = ^(4,57) and X = [4>W[ but X WiW\[i>]- 
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boundary in the context of conformal turbulence while conformal turbulence gen- 
eralizes turbulence to non-inertial regions. In general, the inertial range develops 
regardless the shape of a particular boundary. For example one point functions 
on a disk of radius R can be obtained by finite conformal transformation of Eq.(4) 

[12] 

such that 

This again shows that the minimal dimension operator dominates when R becomes 
large and R <^r. 

Having exemplified the role of one point functions, we note that there is no a 
priori reason to set parameters d^. to zero. It occurs only in the specific case of 
SL{2, C) invariant vacuum or by an unnatural fine tunning. Thus, it is more nat- 
ural to assume that 7^ in our case and use Eq.(lO) for the constant enstrophy 
condition. This together with the Navier-Stokes condition > 2A^ leads us to 
the new set of solutions of conformal turbulence on a upper half plane which are 
listed in the Table 1. In the numerical computation, we have used the following 
notations and rules of 2-d CFT; Mp^q denotes a minimal model with p, q relatively 
prime and p < q, which is characterized by ^(p — — 1) degenerate primary 
operators ^rn,nii ^ n < q, l<m<p)of conformal dimension 

_ {pn - qm)^ - {p - qf 

" ^q ' ^ 

and the fusion rules between two primaries; 

r 8 

^mi,ni X $m2,n2 = XI XI '^(mi,ni)(m2,n2)'^*'-?'' ^^^^ 

i=\m\—m2\+l j=\n-i_—n2\+l 



with r = min(mi-|-m2 — 1, 2p — mi—m2 — 1), s — min(ni-|-n2 — 1, 2g — ni — ^2 — 1), 
and C) ' X/ N are structure coefficients and i, 7 run over odd or even numbers 

(mi,nij(m2,n2) ' 

if it is bounded by odd or even numbers respectively. 
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Finally, a couple of comments related to new solutions are in order. 

It is intriguing to observe that = ^(2,26) occurs repeatedly with increasing 
frequency as we change p and q; 2<p<25, p<g'<500 allof whose energy 
spectrum exponents lie in between —4 and —5. Indeed, there exists a fixed point 
whenp — > 00 and q/p — > 15. In which case ip — '^'(2,26)) ~ ^ ~ ^(1,15)' ^0 ~ 
—49/15; X = ^(2,30)) ^ —49/15, and the energy spectrum exponent is —67/15. 
The physical relevance of this fixed point is unknown. 

Table 1 shows that the energy spectrum exponents can be greater than —3 
whereas in the bulk case they were strictly less than —3. This could be understood 
as a boundary effect to the energy spectrum. Different type of boundaries; e.g. 
different shapes (strip or disk) or periodic boundary conditions can also modify 
the one point functions and the energy spectrum. Details will appear in Ref. [12]. 
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Table 1: Some Solutions for + — = —3 



(p,(l) 




<P 


X 


exponent 









X 


exponent 


(2,33) 


(1,10) 


(1,17) 


(1,16) 


-3.727273 


(3,43) 


(1,12) 


(1,15) 


(1,14) 


-4.837209 


(3,44) 


(2,26) 


(1,15) 


(2,30) 


-4.613636 


(3,46) 


(2,26) 


(1,15) 


(2,30) 


-4.282609 


(3. 52) 


(2. 27) 


(1.17) 


(2. 35) 


-3.307692 


(3. 67) 


(2.31) 


(1.23) 


(2. 15) 


-0.835821 


(3,121) 


(1,12) 


(1,23) 


(1,34) 


-2.000000 


(3,169) 


(1,14) 


(1,27) 


(1,40) 


-2.000000 


(3,196) 


(1,15) 


(1,29) 


(1,43) 


-2.000000 


(3,256) 


(1,17) 


(1,33) 


(1,49) 


-2.000000 


(3,289) 


(1,18) 


(1,35) 


(1,52) 


-2.000000 


(3,361) 


(1,20) 


(1,39) 


(1,58) 


-2.000000 


(3,400) 


(1,21) 


(1,41) 


(1,61) 


-2.000000 


(3, 484) 


(1,23) 


(1,45) 


(1,67) 


-2.000000 


(4,59) 


(2,26) 


(1,15) 


(2,30) 


-4.580508 


(4,61) 


(2,26) 


(1,15) 


(2,30) 


-4.331967 


(4,75) 


(2,28) 


(1,19) 


(2,38) 


-2.590000 


(4,95) 


(3,56) 


(1,23) 


(3,72) 


-0.115789 


(5,71) 


(4,55) 


(1,15) 


(4,57) 


-4.929577 


(5,72) 


(4,55) 


(1,15) 


(4,57) 


-4.777778 


(5,73) 


(2,26) 


(1,15) 


(2,30) 


-4.638356 


(5,74) 


(2,26) 


(1,15) 


(2,30) 


-4.559459 


(5,76) 


(2,26) 


(1,15) 


(2,30) 


-4.360526 


(5,77) 


(2,26) 


(1,15) 


(2,30) 


-4.241558 


(5,91) 


(4, 64) 


(1,19) 


(4, 72) 


-2.890110 


(5, 106) 


(2,30) 


(1,21) 


(2,42) 


-1.371698 


(5,114) 


(4,77) 


(1,23) 


(4,91) 


-0.578947 


(5,124) 


(3,58) 


(1,25) 


(3,74) 


0.424194 
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